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We propose a primitive model of Janus ellipsoids that represent particles with an ellipsoidal core 
and two semi-surfaces coded with dissimilar properties, for example, hydrophobicity and hydrophilic- 
ity, respectively. We investigate the effects of the aspect ratio on the self-assembly morphology and 
dynamical aggregation processes using Monte Carlo simulations. We also discuss certain differences 
between our results and those of earlier results for Janus spheres. In particular, we find that the 
size and structure of the aggregate can be controlled by the aspect ratio. 



INTRODUCTION 

Colloidal particles with anisotropic properties interact 
through an energy that depends not only on their spatial 
separation but also on their relative orientations. This 
is a relatively new field that has been receiving consid- 
erable attention in recent literature pJ-16J. Various site- 
specific techniques such as template-assisted fabrication 
and physical vapor deposition have been developed to 
synthesize patchy colloidal particles of different shapes, 
patterns and functionalities [7-9 . The self-assembly of 
those building blocks into a desired mesoscopic structure 
and function are considered as a bottom-up strategy to 
obtain new bulk materials that have potential applica- 
tions in broad fields including drug delivery, photonic 
crystals, biomaterials and electronics. The effects of 
anisotropy have been classified by Glotzer and Solomon 
using the concept of an anisotropy dimension including 
patchiness, aspect ratio, faceting etc \1\ llQj. 

Patchiness on spheres with different number, size and 
arrangements have been successfully fabricated in col- 
loids experiments and reveal interesting properties and 
crystal structures [SHSl [11]. One particular example is 
the Janus sphere with two dissimilar semi-surfaces, which 
has been extensively studied by experiments and theories 
[3l [T2HT4] . Examples of dissimilar, coded surfaces in- 
clude hydrophobic/hydrophilic, charged/ uncharged, and 
metallic/polymer surfaces. A variety of stable structures 
and unusual phase behaviors have been found under dif- 
ferent chemical conditions of the solution by both exper- 
iments and computer simulations. Without losing the 
generality of the dissimilarity of two surfaces, a primitive 
two-patch Kern-Frenkel model has been used to model 
Janus spheres j2l[T5]. In Monte Carlo simulations, this 
simplified model reproduces the main experimental fea- 
tures including self-assembly morphology and sheds light 
on potential applications in engineering and theoretical 
studies of reentrant phase diagrams [14l Hfi! • 

As suggested by Glotzer and Solomon, it's natural to 
extend the study to explore the role of aspect ratio, in 
which patches are arranged on an anisotropic core such 



as spheroids. The anisotropy dimension is related to the 
aspect ratio and recent studies using ground-state energy 
calculation reveal many interesting structures such as he- 
lix [T7[[lH]. However, how the aspect ratio affects the self- 
assembly morphology is still not clear. In experiments, 
ellipsoidal colloids can be engineered with high monodis- 
persity using techniques such as deforming the spherical 
silica by ion fiuence [I9l [20] , which makes the patched el- 
lipsoidal surface possible if one can combine this with 
the template- assisted fabrication technique. Recently, 
unpublished work shows that Janus football-like ellipsoid 
has been fabricated [21] . Suspensions of Janus ellipsoidal 
particles provide a good candidate to study the effect of 
aspect ratio on self-assembly and new feature of colloidal 
phase transformation. 

In this report, we first propose a theoretical model 
of Janus ellipsoids with hard-core repulsion and quasi- 
square-well attraction. The properties of the model are 
also discussed. In the second part, we use Monte Carlo 
simulations to study the effect of aspect ratio of ellipsoids 
on the self-assembly morphology and dynamical proper- 
ties. In the last section we present a brief conclusion. 

MODEL 

In our simulation study, a primitive model is presented 
for Janus spheroidal particles, with the lengths of the 
principle axes denoted by a 7^ 6 = c. Depending on the 
aspect ratio, defined as e = |, the ellipsoid is charaterized 
as oblate if e < 1 ("M & M") and prolate if e > 1 ("foot- 
ball"). Consider an ellipsoid centered at tq = (a;o,yo,^o) 
whose axial orientations are given by u-^ = {ui^U2^u^)^ 
= {vi^V2^vs)^ W"^ = {wi^W2^ws). The equation for 
this ellipsoid has the explicit form 

[u.(r-ro)]^ [v.(r_i.o)]2 + [w.(r-ro)p ^ 

(1) 

In terms of a matrix representation, Eq. [ijcan be rewrit- 
ten as X-^AX = 0, where X-^ = (x, 1) and A is a 
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4x4 symmetric matrix. The ten independent elements 
of A are listed as follows: 

1) An = g + S + ^ , A12 = ^ + 
Ai3 = + ^ + , = -xoAn 

2) A22 = S + g + #,A23 = ^ 
A24 = -^oA21 - ^oA22 - ^oA23 

3) A33 = ^ + # + ^ , A34 = -xo Ai3 

A44 = -1 - X0A14 - yoA24 - ^oA34. 

Here A is normalized so that a point in the interior of the 
ellipsoid Xq satisfies Xq AXq < and thus det(A) < 0. 
This normalization will be used to determine the spatial 
relation between ellipsoids in the late section. 

Analogous to the patchy spherical model introduced 
by Kern and Frenkel [IS], we choose the ellipsoids to 
interact through a pair-potential that depends on their 
separation and orientation: Uij = [//(r^j, u^, Uj). As 
illustrated for two oblate ellipsoids in Fig. [l] attractive 
patches are coded by red and the orientation is chosen to 
coincide with the principle axis u in the body-fixed frame 
of reference . The orient at ional interaction is defined as: 




where Uq is the well depth, H{x) denotes the Heavi- 
side function, a represents the length of the longer axis 
(max(2a,2b)). Vij is the center- to- center distance be- 
tween ellipsoids: r^.- = |r^j| with r^.- = — r.-. 
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with X = 



e2 + l' 



[22 . aij is introduced as an 



approximation to characterize the spatial relation, such 
that there is no overlapping interaction if r^j > a^j, pro- 
vided that the ellipsoids are represented by a Gaussian 
function: exp(— r •7"-'^ • r) with 7 = a^uu + 6^(vv + ww) 
[22 . This approximation has been widely used to study 
anisotropic particles such as liquid crystals and granular 
material [23]. We note that in our study, this approx- 
imation is only applied to the potential but not to the 
geometric overlapping which will be determined using a 
precise method. Therefore, in our case H{(Jij-\-0.5(j — rij) 
represents a quasi-square- well potential with width 0.5a. 
Specific cases under the condition = Uj are illustrated 
in Fig. [2] for the aspect ratio e = 0.1, 0.5 and 0.9 from 
left to right. Two particles will interact if their red shells 
touch each other. The hard-core repulsion is provided 






FIG. 1: Plot of two interacting ellipsoids labeled as i and j, in 
which attractive and hardcore repulsive surfaces are coded by 
red and blue, respectively. Vij is center-to-center displacement 
pointing from j to i. and Uj is the patchy orientation. Oi 
is the patch angle of ith ellipsoid. 



1 if Ui • Vij < cos (5, Uj • Vij > — cos5 
otherwise 



as shown in Fig. [T] ^ = | corresponds to Janus parti- 
cle. An attractive interaction exists if two red patches 
face each other. The standard square-well potential has 
been used for Janus spheres; however, the determina- 
tion of the accurate spatial relation between ellipsoids is 
computationally time-consuming. We thus introduce a 
quasi-square-well potential defined as 
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if particles overlap 



-UoH{(Jij + 0.5<j — Tij) otherwise 



(3) 



FIG. 2: Plot of quasi- square- well potential under the condi- 
tion Ui = Uj for e = 0.1, 0.5 and 0.9 from left to right. The red 
and blue surfaces represent the attraction range and hardcore 
repulsion associated with each ellipsoid, respectively. View 
angle of each figure is tuned for better visualization. 

by the following geometric relation. Given two ellipsoids 
A: X^AX = and B: X^BX = 0, one introduces the 
characteristic polynomial F{X) = det(A — AB). A and 
B are normalized so that the interiors of A and B satisfy 
X^AX < and X^BX < 0. The roots of the charac- 
teristic equation F{X) =0 have two positive real values 
and the rest characterizes the geometric relation between 
ellipsoids [2Jl [25] such that 

1. A and B are separate if and only if F{X) = has two 
distinct negative roots; 

2. A and B touch each other externally if and only if 
F(A) = has a negative double root. 

3. Otherwise, A and B overlap. 

Sturm sequence methods are applied to numerically de- 
cide if two roots are distinguishable. The primitive model 
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we propose recovers the Kern-Frenkel Janus sphere model 
when e = 1, since under this condition, aij = a and con- 
sequently, the attraction is simplified to H{l.ba — Vij). 
The introduced quasi-square- well potential has an advan- 
tage that there is no ambiguity when defining the con- 
nectivity of aggregates during the self-assembly process, 
and could be easily generalized to more realistic models. 

Our model is different from the Kern-Frenkel model 
when e 7^ 1 since then the magnitude /7 is a func- 
tion of both the separation and orientation due to the 
anisotropic ellipsoidal core. The aspect ratio e affects 
the shape of the interacting potentials. To characterize 
the effective interaction between particles, we calculate 
the second viral coefficient B2'. 

B2 = /[I - e"^""^^] dridr2duidu2 (5) 

The results of Monte Carlo integrations of B2 / ^2 * as a 
function of temperature for different e are shown in Fig. 
[3| where without specification, e = 0.1,0.5 and 0.9 are 
denoted by green down triangles, red triangles and blue 
circles, respectively. Here ^2 ^ = ^Tra^ stands for the sec- 
ond viral coefficient for hard spheres with diameter equal 
to <j. The numerical error is less than 1%. The solid curve 
represents the theoretical prediction for a Janus sphere: 
B2/B^' = 1- 1(6^ - l)(e^^o - 1) with the interaction 
range 6 = 1.5a in the study [15]. As e increases to 1, B2 
approaches the value for the Janus sphere from above. 
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FIG. 3: Plot of vs. kBT/Uo for the aspect ratio e = 
0.1 (green down triangle), 0.5 (red triangle) and 0.9 (blue 
circle) from top to bottom. The solid curve is the theoretical 
prediction for a Janus sphere, i.e. e = 1. 



ary conditions. The particle aspect ratios range from 
0.1 to 0.9 with number density p = 0.037; the attractive 
energy is set as /3Uo = 3 with /3 = where ks is 

the Boltzmann constant. The system is initialized as a 
randomly-distributed noninter acting gas of monomers. A 
random translation followed by a random rotation is car- 
ried out for each monomer. In particular, the rotation is 
performed using the method of quaternion parameters. 
More than ten independent runs for each aspect ratio 
have been carried out up to 5 x 10^ Monte Carlo steps 
(MCS) and an ensemble average is taken by averaging 
over all runs. We monitor the evolution of —E/Uq (shown 
in Fig. [4|, the negative average potential scaled by the 
attraction strength, which characterizes the number of 
interacting neighbors of each particle. Since the initial 
system has no interaction, —E/Uq starts from 0. It grows 
quickly and the dynamics slows down while approaching 
equilibrium. The system for e = 0.9 has reached equilib- 
rium; however, for e = 0.1 and 0.5 the system is only close 
to the equilibrium state after 5 x 10^ MCS. This differ- 
ent growth tendency for the different aspect ratios is due 
to the distinct aggregation mechanism that dominates at 
different times. As we will show in the later section, the 
aggregation at the early stage is dominated by monomer 
diffusion and interactions with small, formed oligomers. 
Depending on the aspect ratio, the system is composed 
of monomers, small oligomers, micelles and vesicles. Af- 
terwards, the aggregation dynamics is mainly the diffu- 
sion and collision of those small clusters that have much 
smaller diffusion constants than monomers. The value of 
e affects the dynamics, as shown in Fig. |4j such that the 
aggregation process is relatively faster but reaches a less 
stable structure when e is larger. 
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SIMULATION RESULTS 

Standard Monte Carlo (MC) simulation in the NVT 
ensemble has been applied to study the self-assembly 
pathway of these interacting Janus particles. We investi- 
gate a system in a 30 x 30 x 30 box with periodic bound- 



FIG. 4: Plot of -E/Uo vs. MCS. The curves from top to 
bottom correspond to the aspect ratio e = 0.1 (green), 0.5 
(red), and 0.9 (blue) 

To illustrate the difference in aggregation morphology, 
we take snapshots of the cluster growth as shown in Fig. 
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FIG. 5: Plot of cluster formation for e = 0.9, 0.5 and 0.1 from left to right. The number next to each panel is the cluster size. 
The size and view angle of each figure have been tuned for better visualization. 



m The cluster is defined such that two monomers are 
connected if they interact; there is no ambiguity with re- 
spect to this interaction in our model. From left to right 
in Fig. |5) we show typical structures for e = 0.1, 0.5, 
and 0.9 that are collected from all simulations. Small 
oligomers such as trimer, tetramer, pentamer, hexamer, 
and heptamer are similar and monomers form usual poly- 
gons, except that for e = 0.1 there is more space enclosed 
by monomers and the heptamer has a different structure: 
a rough double-layer. When the aspect ratio is less than 
0.5, a different feature shows up. Instead of building up 
polygons, two ellipsoids with the same orientation pack 
into two layers. This is due to the quasi-square- well at- 
traction in our model of the attraction range 0.5<j. When 
e < 0.5, two ellipsoids have the possibility to interact 
even when they are separated by the third one. As the 
oligomer grows, it starts forming a micelle (single layer) 
and a vesicle (double layer) as shown for clusters 13 and 
39 (e = 0.9), and 8 and 16 (e = 0.5). For e = 0.1, the 
relevant clusters (8 and 17) display two and three layers, 
respectively. Those structures have been found in the 
case of Janus spheres and have an effect on breaking the 
thermal correlations between particles. They thus affect 
the system phase behavior. Eventually for e = 0.9, clus- 
ters with different structure are formed, including two 
small oligomers joining together (46), as well as more 
complicated chains (72). For e = 0.5, it's possible to 



form a triple-layer compact cluster (51) and a dumbbell 
cluster (62). For e = 0.1, due to the relation between the 
attractive range and aspect ratio, multiple-layer struc- 
tures have been observed in the simulations (30, 52). 

As is well known, Monte Carlo simulation doesn't pro- 
vide an accurate description of the kinetics of the sys- 
tem; however, it still reveals some useful features. We 
investigate the time evolution of number of clusters as 
illustrated in Fig. [6j The number of clusters for larger 
e drops faster until reaching about 10^ MCS. After that, 
the tendency reverses. The observation is consistent with 
the energy evolution and B2. At the early stage, the dy- 
namics is dominated by monomer motion and the sys- 
tem at e = 0.9 has a relatively stronger interaction which 
leads to faster cluster formation. Then the dynamics is 
governed by cluster-cluster interaction so that for larger 
aspect ratio the number of clusters decreases more slowly 
than for smaller aspect ratio. 

Note that the cluster defined by energy interaction 
might in principle be inconsistent with the experimen- 
tal observations, which could implicitly use the distance 
definition without considering orientations. The inset in 
Fig. [6] shows that the difference between the energy def- 
inition (upper blue) and distance definition (lower red) 
for e = 0.9 is smaller than the statistical error, which 
indicates that the energy-defined cluster is a good ap- 
proximation as compared with experiments. 
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FIG. 6: Plot of number of cluster vs. MCS for e = 0.9 (blue 
circle), 0.5 (red up triangle) and 0.1 (green down triangle). 
Error bars come from the statistical variance of independent 
runs. The inset shows the evolution of the number of clusters 
defined through energy (blue) and distance (red) for e = 0.9. 

We have calculated the distribution of cluster size, 
which illustrates the effect of aspect ratio on the system 
approaching equilibrium, as shown in Fig. [7[ For larger 
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FIG. 7: Plot of distribution of cluster size for e = 0.1 (green 
down triangle), 0.5 (red triangle) and 0.9 (blue circle). Insets 
from top to bottom are configurations of the largest cluster 
found in the simulations for e = 0.1, 0.5 and 0.9, respectively. 
The number associated with each configuration is the cluster 
size. 

e, the distribution has a broader range with a lower peak, 
which is consistent with the morphology shown in Fig. |5j 
Due to the complex structure formation such as a chain, 
6 = 0.9 has extended configurations with larger size, in 
which the simulations show the largest cluster (a size 137 
as shown in the inset of Fig. [t]). For e = 0.1, the distribu- 
tion is narrow and has a higher peak, which indicates that 
the clusters are more uniform. The largest cluster with 



size 55 has a similar shape as the oligomer (52) shown 
in Fig. [5] the structure of this cluster is like a blob with 
several layers. The structure is more stable in terms of 
its energy and prevents further aggregation to form more 
complex forms. The case of an aspect ratio 0.5 (74) is 
intermediate between these two extremes of 0.1 and 0.9 
and the resulting cluster reveals two vesicles forming to- 
gether through a bridge-like structure. Next, in order to 
investigate in detail how ellipsoids organize in the cluster, 
we consider the correlation between patch orientations: 
Ui • Uj of two bonded ellipsoids (for which there exists a 
patchy attraction). The distribution P{ui • Uj) is illus- 
trated in Fig. jS] 




VLi'Ui 

FIG. 8: Plot of distribution of • uj between two bonded 
ellipsoids for e = 0.1 (green), 0.5 (red) and 0.9 (blue). Typical 
configurations corresponding to each peak are shown. 

For e = 0.9, P(u^, Uj) develops two peaks; typical con- 
figurations corresponding to this correlation are shown in 
the figure. These two peaks correspond to two ellipsoids 
facing toward (u^ • uj ~ — 1) and opposite to (u^ • uj ~ 1) 
each other; similar behavior has been found for Janus 
spheres [12]. As e decreases, one peak ( • Uj ~ —1) 
vanishes, since ellipsoids self-assemble into a micelle-like 
structure (for example, panels 8 and 16 for e = 0.5 and 
panels 8 and 17 for e = 0.1 in Fig. [5]). Due to the geomet- 
ric constraint, there is empty space inside the ellipsoids 
forming the cluster. This property could be useful for 
the encapsulation of different types of particles (e.g. as 
in drug delivery). The other peak shifts towards zero; 
thus, the curvature is increasing so that the structure is 
more compact and energetically favorable. This structure 
prevents further aggregation, as the thermal fluctuation 
is suppressed, which prevents long range correlation. 

More details could be explored from the snapshots of 
configurations. As illustrated in Fig. |9) the system for 
e = 0.9 forms large clusters at 10^ MCS and is composed 
of a large number of monomers, micelles, vesicles and 
large complex clusters at 5 x 10^ MCS. 

For e = 0.5, the configuration shows a somewhat dif- 
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FIG. 9: Snapshots of the system for e = 0.9, 0.5 and 0.1 from left to right. The upper and lower panels correspond to t = 10^ 
and 5 X 10^ MCS. 



ferent feature (middle panels of Fig. |9| . The system 
forms micelles and vesicles at 10^ MCS and finally in- 
cludes monomers, large numbers of micelles and vesicles, 
and their combined aggregates. For e = 0.1, the configu- 
ration show that clusters forms slowly and only includes 
small oligomers and a large number of micelles and vesi- 
cles at 5 X 10^. 



CONCLUSION 

In the article, we propose a primitive model to study 
the self-assembly of Janus ellipsoids. The interactions 
between ellipsoids include a hard-core repulsion and a 
quasi-square- well attraction, where the latter exists when 
the patchy surfaces of two interacting ellipsoids orient 
in designated directions. The anisotropy in our model 
comes from two aspects: the patch interaction and the 
anisotropic core, which are controlled by a patchy angle 
S and an aspect ratio e, respectively. We particularly fo- 
cus on the Janus ellipsoids in which S = ^ and e ranges 
from 0.1 to 0.9, and address the effects of aspect ratio on 
the self-assembly morphology and dynamical properties. 
Our model could be easily extended to consider more 
realistic systems. For example, the typical interaction 
between colloids is in the range (0.05cr ~ 0.2cr), which 



however is computationally expensive to simulate. Janus 
ellipsoids with the interacting range comparable to 0.5<j 
can be realized by nanoparticles using the technique such 
as induced phase separation [26 . Our results show that 
for larger aspect ratio, the self-assembly process is rel- 
atively faster and the morphology is more complicated, 
including chain structures as well as micelles and vesicles. 
The structures for smaller aspect ratio are more uniform 
and multiple-layer vesicles dominate. In our Monte Carlo 
simulation, for e = 0.1 and 0.5, the system didn't reach 
equilibrium since when multiple-layer vesicles form, the 
simulation is extremely time-consuming. However, the 
main features do not change, since the systems are very 
close to equilibrium. In our study, the phase diagram 
is not yet known, but our initial choice of parameters 
(such as density and temperature) puts us in the gas 
regime. Our simple model raises many potential avenues 
for investigation; for example, the possible extension of 
52-scaling in the calculation of the phase diagram. This 
B2 scaling has been successfully applied, for example, for 
patchy spheres [27 . In addition, it is important see how 
our results depend on the range of the interaction, e.g. 
0.2(7, as is the case for colloidal interactions. We also find 
for small aspect ratio that ellipsoids tend to form vesi- 
cles, which suggests a potential application for particle 
encapsulation. Perhaps the most important result of the 
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study for materials engineering is the fact that the size 
and structure of the aggregates can be controhed by the 
aspect ratio, which should be an interesting result from 
a design viewpoint. 
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